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A recent simplified transfer matrix solution of the two-dimensional Ising model on a square lattice 
with periodic boundary conditions is generalized to periodic-antiperiodic, antiperiodic-periodic and 
antiperiodic-antiperiodic boundary conditions. It is suggested to employ linear combinations of the 
resulting partition functions to investigate finite-size scaling. An exact relation of such a combination 
to the partition function corresponding to Brascamp-Kunz boundary conditions is found. 
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Recently, a solution of the two-dimensional Ising 
Model on a square lattice with periodic boundary- 
conditions was presented [jjj, which simplifies Kauf- 
man's version || of Onsagcr's solution || without 
requiring sophisticated mathematics. Here we point 
out how the method easily generalizes also to the 
cases of periodic-antiperiodic, antiperiodic-periodic and 
antiperiodic-antiperiodic boundary conditions, which 
were recently solved using a Grassmann path integral 
approach p|. The motivation for varying the bound- 
ary conditions originates in the investigation of finite- 
size scaling (FSS) [H [| Q]- The analytic properties of 
FSS are of interest when trying to describe critical prop- 
erties of models whose exact solution is unknown. For 
the analysis of FSS in solvable models it is very help- 
ful if the exact partition function can be written as a 
product rather than as a sum ^ Q as is the case with 
Brascamp-Kunz (BK) boundary conditions ||. Here we 
exhibit linear combinations of partition functions with 
varying boundary conditions that have this property and 
may be used for such an analysis. For one of these linear 
combinations we find an exact relation to the partition 
function for BK boundary conditions. 

We follow the notation in Q|. Define the energy E of 
the two-dimensional Ising model on a square lattice with 
to x n sites by 

- f3E = a}^ I <7 a s mv si v + s^s^+i^ J 

m / n— 1 \ 

, (1) 

M =l V v=l / 

where we assume a, b > as in Q and where the spin 
variables s^ v can assume the values ±1. (cr^erf,) = 
(1,1) corresponds to periodic boundary conditions in 



both directions (pp), as investigated in [[j], ||, 
(ca,er;>) = (1,-1), (—1,1) and (—1,-1) correspond to 
periodic-antiperiodic (pa), antiperiodic-periodic (ap) and 
antiperiodic-antiperiodic (aa) boundary conditions, re- 
spectively. The partition function is defined by 



(2) 



s u „=±l 



where af3 may be pp, pa, ap or aa. After defining 2™ x 2™ 
matrices X„ (and analogously Y„ and Z„) by 

(3) 



where a x is a Pauli matrix and 1 is the 2x2 unit matrix, 
we may write the partition function as 

z ^' n \a, b) = [2sinh(2a)] mn/2 Tr(QV m ) (4) 

where Q = 1 (now the 2™ x 2™ unit matrix) and Q = 
Ux = X1X2 ■ ■ ■ X„ for periodic and antiperiodic boundary 
conditions in the a direction, respectively, and where 



v = v Q/2 v b v a/2 



(5) 



with 



where a is defined by sinh 2a sinh 2a = 1 as in (!]] , and 

n-l 

V b = exp(o- b oZ„Zi) Y[ cxp (bZ u Z v+ i) (7) 



with CTb = 1 and at, = — 1 for periodic and antiperi- 
odic boundary conditions in the b direction, respectively. 
With the same set of matrices 
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J a/3 



P±J 



1 



a/3, 



(1±U X ) 



(8) 



2 



as in Q, where 



Ja/3 

r2i/-i = Xi • • • Xj,_!Z„, 

T2i/ = Xi ■ • ■ X^-iY^, 
the partition function may be rewritten as 

^S' H) M) = [2sinh(2a)] m "/ 2 Tr(Q Q V™) 

with 

Q { p } = P+±P_, 

y = v+v^ =v^v 



/3 v /3 ' 



(9) 
(10) 
(11) 



(12) 

(13) 
(14) 
(15) 



half of the terms in the partition function, see (|32j) and 
© below. 
Define 

X£' n) (a,b) = E ex p(f EW^J, (20) 
&™> n Ha,b) = ^ex P (^f](±) 72 _ 1 ), (21) 



E&' w) (a,6) = ^exp(^(±) 72 



2i/-l 



(22) 



(23) 



where the first index in S xy and the index under the 
summation sign refer to all combinations in the sum in 
the exponent that have an even/odd number of minus 
signs and the second index of S xy refers to even/odd 
indices of the 7 fc. Define also 



y %2 = II^^^-i) 



n-l 



u=l 



(16) 



V± = exp( T 26J± 2 „) J] exp(26J± +12 J, (17) 



V± = exp(±26J± 2 „) J] exp(26J± +li2 J. (18) 



From here on it is straightforward to follow the argu- 
ments given in Q to compute the partition function in 
terms of the 7& defined by 

cosh7fc = cosh 2a cosh 2b — cos — sinh2asinh26 (19) 

n 

with the understanding that jk > for 1 < k < n and 
with 70 = 2 (a — b). It turns out that antiperiodic bound- 
ary conditions in the 6 direction lead to an exchange of 
7^ with even and odd indices, while antiperiodic bound- 
ary conditions in the a direction lead to a sign change of 



C^ n) (a,b) = [2sinh(2a)] m ™/ 2 jj [2cosh(y 

fc=i 



7 2fc-l 



(24) 



n 

Si m ^(a,b) = [2sinh(2a)] mn / 2 [] [2 sinh 



fe=i 



(25) 



n 

C^ m ' n) {a,b) = [2sinh(2a)]" m / 2 ]J[ [2 cosh (— 



72fc-2 



fe=l 



(26) 



S< m ' n )(a,&) = [2sinh(2a)] m ™/ 2 J| [2sinh(^ 

fc=i 



and note that 



72fe-2 



(27) 



Sf ro '")(a,6) = S { m ' n) (a,b)C^ n \a,b), (28) 
S( 2ro '")(a,6) = 5f l '" ) (a,6)Ci m ^ ) (a,6). (29) 

Then, omitting upper indices (m, n) and arguments (a, 6) 
for brevity, we have for the various cases 



Z pp = [2sinh(2a)] mn/2 (E co + S oc ) = -( 

Z pa = [2sinh(2a)]"" l / 2 (E cc + S 00 ) = i( 

Z ap - [2sinh(2a)] mn / 2 (E co -S oc ) = i( 

Z aa = [2sinh(2a)] mn / 2 (E cc - E DO ) = i( 



Co + S'o + Ce-S'e), (30) 

Co-S'o + Ce + S'e), (31) 

C + S -C c + S c ), (32) 

Co + S'o + Ce + S'e). (33) 
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The result for Z pp is the one found in M and, observing 
the different sign convention for 70, in p|. 
In §, K FSS 

was investigated with BK boundary 
conditions |8j so that the partition function has a prod- 
uct structure which facilitates the subsequent analysis of 
the approach to the critical point in the thermodynamic 
limit. The partition functions above do not have such a 
product structure in contrast to their linear combinations 
C , iS OI C c and S c . Except for S e , these combinations are 
positive for any temperature and may be viewed as parti- 
tion functions belonging to models with possibly nonlocal 
boundary conditions. Even though these boundary con- 
ditions may in general not be particularly physical, the 
quantities C , 5*0 and C e may be used as a mathemati- 
cal tool to investigate the approach to criticality in the 



same way as the partition function of the model with BK 
boundary conditions. 

S e does not qualify for such an approach. Since the 
pp boundary conditions cause a lower ground state en- 
ergy than the other three boundary conditions considered 
here, S c becomes negative for sufficiently low tempera- 
tures since Z pp enters with a minus sign. S e has another 
interesting property, though: Due to the sign change of 70 
at the critical temperature of the thermodynamic limit, 
(i.e., where a = b), S c changes its sign at the critical 
temperature for any size mxnof the square lattice. 

The quantities (p4|)-(p7|) are the same as in even 
though they have a different appearance. In our notation, 
the expressions in § read 



1 " ' (2 + l)ir (2 + l)ir 1 1 / 2 
cosh(2a) cosh(26) — cos sinh(2a) — cos sinh(2&) . 



C^ m ' n \a,b) = 2 mn Y[ II 

p=0 q=0 
m— 1 n — 1 

S { m > n) (a,b) = 2 mn \\ Y[ 

p=0 q=0 
m-ln-1 r 

C^ m ' n) (a,b) = 2 mn Y[ ]J 

p=0 q=0 

S { m > n) (a,b) = sgn(l - sinh(2a) sinh(26)) x 

ml n—1 r 

2 mn y[ n 



cosh(2o) cosh(26) — cos sinh(2a) — cos - - — — sinh(2&) 
777 n 

cosh(2a) cosh(26) — cos — — — - sinh(2a) — cos sinh(2&) 

m n 



p=0 q=0 



2j)7T 2q7T 

cosh(2a) cosh(26) — cos sinh(2a) — cos sinh(2fr) 

m n 



1/2 



1/2 



1/2 



(34) 
(35) 
(36) 

(37) 



While these expressions involve a double product instead 
of just a single product as in (24)-(p71), they have the 
advantage that C and S e are symmetric under the ex- 
change (a, to) <-> (6, n) and that the only asymmetry of 
C e and S under this exchange originates in the boundary 
conditions and not in the mathematical treatment. 

While Wu et al. j§ use the Grassmannian integral ap- 
proach to map the problem onto a free fermion system, 
the method employed in this work is close in spirit to 
Kaufman's solution || , but avoiding some complications 
encountered there, as detailed in While the Grass- 
mannian method may be viewed as more elegant, our 
method is accessible at a less sophisticated mathemati- 
cal level. It is also interesting to note the very different 
appearance of the solutions (PI)-© and It 



would be interesting to see if our method can be gen- 
eralized to cases that so far have been solved only us- 
ing the Grassmannian integral approach, such as various 
boundary conditions on plane-triangular and honeycomb 
lattices ||, or other. 

Brascamp and Kunz Q gave a solution for a two- 
dimensional Ising model on a square lattice with (771+2) x 
77 sites with periodic boundary conditions in the b direc- 
tion and fixed boundary conditions in the a direction. If 
k enumerates the sites in the a direction (0 < k < to+1), 
then the spins at k = are fixed at +1, while the spins 
at k = 777+1 are fixed to alternate between —1 and +1, 
so 77 must be even. The resulting partition function was 
determined to be Hi 



n/2 r 



p=l q=l 



J)7T (.^Q — lj7T 

cosh(2a) cosh(26) — cos sinh(2a) — cos sinh(26) 

777+I 77 



(38) 



r 



Using ( pig ) and (p4|)-(p6D, we can establish the connection 



= Ct»{b,a)Z^ n \a,b)\ (39) 



4 



where the S Q may be factorized according to (|2S|) and 
ci n,2 \b 2 a) most easily evaluated using (p6|). 

Eq. (||) generalizes Eq. (10) in @ to the case a ^ b. It 
would be interesting to understand ( |39| ) without having 
to resort to the solutions of the models. Then it may 
also be possible to relate C , S , C c and S e to yet other 
physical boundary conditions. 

Investigation of FSS for C a , S a and C c is beyond the 
scope of this paper. Let us however state one interesting 
result that can be read off immediately from ([34|). If 
a = b and we define z = sinh2a, the zeros z pq of C Q in 
the complex plane are given by z pq — exp(ia P9 ) with 



arccos 



1 / (2p+l)7T 

- cos 

2 V m 



cos 



(2q+ 1)tt 



(40) 



and the zero nearest to the real axes is found for p = q 
0. For m — n, this zero scales as 



7T 

«00 = - 

n 



(41) 



without any corrections as opposed to the case of BK 
boundary conditions [HQ]. It remains to be investigated 
if simplifications may also be achieved for the FSS of 
other quantities and if simplifying combinations of parti- 
tion functions can also be defined for models whose exact 
solution is unavailable, e.g. by the use of symmetry ar- 
guments. 
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